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ABSTRACT — The design, analysis, and demonstration of a digital-correlation microwave 
polarimeter for use in earth remote sensing is presented. We begin with an analysis of three- 
level digital correlation and develop the correlator transfer function and radiometric sensitivity. 
A fifth-order polynomial regression is derived for inverting the digital correlation coefficient into 
the analog statistic. In addition, the effects of quantizer threshold asymmetry and hysteresis are 
discussed. A two-look unpolarized calibration scheme is developed for identifying correlation 
offsets. The developed theory and calibration method are verified using a 10.7 GHz and a 37.0 GHz 
polarimeter. The polarimeters are based upon 1 -GS/s three-level digital correlators and measure 
the first three Stokes parameters. Through experiment, the radiometric sensitivity is shown to 
approach the theoretical as derived earlier in the paper and the two-look unpolarized calibration 
method is successfully compared with results using a polarimetric scheme. Finally, sample data 
from an aircraft experiment demonstrates that the polarimeter is highly-useful for ocean wind- 


vector measurement. 



I Introduction 


Recent advances in the interpretation of polarimetric microwave thermal emission from the Earth’s 
oceans and atmosphere have prompted the study of new retrieval techniques for near-surface ocean 
wind vectors and mesospheric temperature profiles [1], These techniques are facilitated by a more 
complete characterization of the polarization characteristics of the upwelling radiation field than 
obtainable using conventional single- or dual-polarization radiometers. As an example of these 

techniques, polanmetnc measurements have been shown to greatly facilitate the retrieval of ocean 
surface wind direction [2,3]. 

The quantity used to fully describe the second-order statistics of the quasi-monochromatic 
radiation field at a point in space is the Stokes vector. The elements of the modified Stokes vector, 
in units of brightness temperature (Kelvins), are directly related to the following ensemble averages 
of the incident transverse electric field components [4]: 
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where A is the wavelength, rj is the impedance of the medium, and k is Boltzmann’s constant. Here, 
E a is the phasor amplitude for polarization a (= v or h) per unit solid angle and bandwidth. The 
first two parameters, T v and T h , are the intensity in the vertical and horizontal polarizations and 
their sum is the total radiation intensity. The remaining two parameters contain information about 
the polarization characteristics of the radiation field. Specifically, T v indicates the degree and 
sense of linear polarization and T v of circular polarization. Partially polarized thermal radiation is 
specified by nonzero T v or T v and by T v ± T h . 

The parameters T v and T h can be measured using standard linearly-polarized total power 
radiometers [5], Detection of the third and fourth Stokes parameters, however, requires two ad- 
ditional measurements to effectively perform the correlations in (1). The various types of polari- 
metnc radiometers fall into two basic categories: adding polarimeters (AP) and direct correlating 
polarimeters (DCP). The adding polarimeter uses measurements of the brightness temperature of 
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at least two additional polarization states e.g., 45° slant-linearly polarized (T 45 ») and either lefit- 
or right-hand circularly polanzed (7] or T r ). From the four measured brightness temperatures and 
using the Stokes parameter rotational transformation [6], the third and fourth Stokes parameters 
can be determined, e.g.: 


Tu — 2T 45 o — T v — T h (2) 

T v = 2T r -T v -T h (3 ) 

The direct correlating polarimeter estimates T v and T v by cross-correlating the instanta- 
neous voltage signals of the vertical and horizontal channels. The actual correlation can be per- 
formed by either analog or digital multiplying circuitry. If the time-varying voltages v v (t) and 
v h (t) are assumed to be stationary and ergodic [7], then the covariance estimate Ry h is: 



f 


v v {t)v h (t)dt 


( 4 ) 


where r is the integration time. Since the IF voltages are related to the incident field quantities by 
the receiving antenna’s effective area and the receiver’s signal transfer characteristics, measuring 
Rvh is equivalent to measuring Tu: 


Tu = 2 py/f Vity ,f, 


h,sys 


( 5 ) 


where p s v d h * s ^e correlation coefficient and T a sya are the system temperatures of the ra- 
diometers (assumed to be of the total power type) for polarizations a = v and h. A dual-channel 
superheterodyne receiver with a coherent local oscillator (LO) may be required to downconvert 
the RF band of interest to accommodate the bandwidth and/or operating frequency of the analog 
multiplier or the digital correlator A/D converters. If the signals are downconverted in-phase, T v is 
estimated; however, if the receiver is single-sideband and the signals are downconverted in phase- 

quadrature, then T v is estimated. Direct correlation without down conversion is also possible with 
the appropriate hardware. 

Several mechanisms can contribute to calibration errors in (4) and (5). Antenna cross- 
polarization mixing and receiver imbalances can cause mixing of all four Stokes parameters, the 
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amount of which must be known. One method for comprehensive calibration of the first three mod- 
ified Stokes parameters uses a rotating polarized calibration standard [8]. The polarized standard 
presents to the receiver a strongly polarized but precisely determined radiation field and allows 
complete determination of the gains and offsets for the first three Stokes parameters. Calibration 
of the fourth Stokes parameter channel can be accomplished by insertion of an appropriate 90° shift 
in the RF path using, e.g., a quarter wave plate [9]. Use of the polarized standard in space, however, 
requires additional hardware beyond the conventional ambient and cold blackbody standards that 
are commonly used. 

In the implementation of (4) and (5) it is desirable to design a system that requires a minimal 
amount of calibration hardware. While an analog correlator can be used to determine Ty or Ty 9 its 
response will generally require the in-situ identification of relatively large leakage gains g Uv and 
guh from T v and T h , viz.: 


v u — 9 uvT v + guhTh + guuTu + guvTy + ojj ( 6 ) 

as well as the offset term ou . While leakage gains can be minimized by proper tuning and balanc- 
ing, elimination of long term drift in detection and video components - their root cause - can be 
prohibitively expensive. 

A solution to the above problem above of precise measurement of either Tjj or Ty can be 
found through digital correlation. Here the radio-frequency (RF) or intermediate-frequency (IF) 
signals are sampled at the Nyquist rate, the digital samples cross-correlated using fast multiplica- 
tion circuitry, and the correlation integral (4) performed via digital accumulation. Provided that 
the digitized signal contains no DC component and the A/D conversion is linear and unbiased, the 
correlation coefficient p can be obtained without leakage or offset. A further advantage of using a 
digital correlator with more than one bit (or two levels) of discretization is that in-situ calibration 
can be performed using only conventional ambient and cold unpolarized views, for example, an 
ambient blackbody target and cold space. 

Digital correlation radiometry was first discussed by Weinreb [10] for use in autocorrela- 
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tion spectrometers for radio astronomy. As shown by Weinreb in autocorrelation spectroscopy of 
Gaussian signals, only a single bit of quantization (i.e., two-level A/D conversion) is required to 
achieve ~64% of the detection sensitivity of a perfect analog correlator system. As fast digital 
logic became more widely available the single-bit systems were replaced with three-level (reduced 
2-bit or 1.6-bit) systems. The three-level correlator can obtain up to 81% of the detection sensi- 
tivity of the analog system [11], The availability of discrete high-speed digital logic has allowed 
spectrometers to operate over wider bandwidths, and both single-bit and two-bit correlators have 
now been implemented at clock-rates as high as 2 GS/s (e.g., [12-16]). 

In this paper we present the first demonstration of a polarimetric digital radiometer for earth 
remote sensing. Because of the relatively wide bandwidths required for earth remote sensing appli- 
cations (typically tens to thousands of MHz), the digital correlator has not been considered for use 
in space until recently. With the advent of high-speed radiation-hardened digital logic, bandwidths 
of hundreds of MHz have now become detectable using sensitive, short integration-time, digital 
correlating polarimeters. The block diagram of a generic digital correlating polarimeter is shown 
in Figure 1. The major components are a (i) dual-polarized antenna, (ii) a superheterodyne phase- 
coherent dual-channel receiver, and (iii) a three-level digital cross-correlator). We begin in Section 
II with a description of digital correlation radiometry and discuss in Section III an investigation 
of systematic errors along with design implications. A novel two-look calibration of the correla- 
tion channel is discussed in Section IV. The results obtained from an airborne 1.6-bit X-band and 
K a -band 1 GS/s digital polarimetric radiometer are reported in Section V, thus demonstrating the 
ability to properly calibrate and successfully operate a digital polarimetric radiometer. 

II Digital Correlation Polarimeter 

The digital direct-correlating polarimeter has as its main distinguishing component a zero-lag dig- 
ital cross-correlator. The digital correlator is composed of three sections: a pair of A/D converters, 
a digital multiplier, and accumulators. To fully understand the correlator operation, it is necessary 
to examine the relationship between the input signal statistics and the accumulator outputs. These 
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relationships also provide a measure of the digital correlator sensitivity. 

A Mean Statistics 


The input signals to a correlator, v a (t) and v b (t), are modeled as jointly-Gaussian stationary random 
processes with root mean square (RMS) voltages o Va and a v and correlation coefficient o = — 

If the processes are sampled with period T at or below their Nyquist rate, then the sample sequences 

consist of independently and identically distributed pairs with the following joint Gaussian proba- 
bility density function (pdf): 


f{v a ,v b ;p ) = 


1 

flvn 
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The three-level quantization performed on the input signals by the A/D converter can be modeled 
by a nonlinear transfer function: 


I 1 if v > v tha 

-1 iftK-t^ (8) 

0 otherwise 

where the quantities ±v tha are the threshold levels of the A/D converters (see also Figure 2), with 
the subscript a denoting either channel a or b. For typical CMOS or ECL logic, v tha ^ 0.05 to 
0.50 volts, therefore, the necessary microwave signal power can range from -12 to +8 dBm in a 
standard 50 $1 system. 


The outputs of the quantizers form a new pair ofjointly-random processes denoted h[y a [nT )) 
and h(v b (nT )) where sample n is taken at time nT. The second-order statistics of these sampled 
and quantized processes are the digital variances and covariance and are nonlinearly related to 
the first three Stokes parameters. For a measured sequence of N samples, the estimated digital 
variances and covariance, denoted §1 and f ab , are: 

1 N 

Q = jyX^MnT)) 2 (9) 


1 

fab = 

n — 1 


( 10 ) 
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These three statistical parameters are measured by accumulation of the outputs of a simple digital 
circuit such as shown in Figure 3. 

The statistics of sjj, s b , and f ab and their relationship to T a , T b , and T v are obtained by 
integrating the right-hand sides of (9) and (10) against the pdf (7). The expected value of the 
digital variance is 


<£> = 2 [1 - •(*„)] 


( 11 ) 


where 9 a = v tha /cr Va and 


*(*) = 4= r 

V2ttJ-oo 


( 12 ) 


is the normal cumulative distribution function. Figure 4 is a plot showing the relationship between 
the digital variance and RMS input voltage at a fixed threshold voltage. As will be shown in section 
B, for maximum sensitivity in T v the value of 9 a should be close to 0.61. Inverting (11) yields a 
simple estimate for the signal standard deviation for a measured digital signal variance: 


k = 0)- 1 



(13) 


or, in terms of antenna brightness temperatures: 
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where Ro is the system impedance, D is the bandwidth, G a is the system gain, and T RECya is 
the receiver noise temperature. In general, the parameters ( B tt Ga ) and T RECa are slowly time 
varying and represent system gains and offsets that must be identified via periodic calibration. 

The relationship between the input correlation coefficient p and the expected value of the 
digital covariance r ab is similarly straightforward and can be obtained by integrating the right- 
hand side of (10) against the joint pdf over two dimensions. The problem can be reduced, however, 
to an integration over one-dimension using Price’s theorem [17, 18]. Price’s theorem relates the 
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covariance of the input signals to the digital correlation coefficient: 
drgb _ / dh(v a ) dh(v b ) \ 

dRv aVb \ dv a dv b J 

= <[£K + v tha ) + 6(v* - v th a )] + v th b ) + ${v b - w«fc»)]> 

= f ( v thll , v thb ;/>) + / (««* , -v thb ; p) + /(-uth 0 , ; p) + f{-v tha , -v t h b ; p) ^ 1 

= /K iA, ^„ 6 6> 6 ; p) + f(-a v J a , o Vh 0 b \ p ) + 

/(<7„ o 6> 0 > —Ov b Qb\ P) + f(-CT Va e a, -0 Vi 9 b ] p) 

where J is the Dirac impulse function. The input covariance can be related to the input correlation 
coefficient using the chain rule of differentiation: 

d'f'ab 9t a i dR Va Vb dr a b ,, 

“q on O "qd ( 1 h) 

up uRy aV{} op dR VaVb 

Thus, the digital correlation coefficient is a one-dimensional integral of the pdf over p: 
fab = &v a @v b f [f{<?v a O a , &v b 9b > P ) / (°i ; a $o> <fv b @b: P ) T 

Jo 

f{-o Va 0a, cr v „9 b , p') + f(~a Va 9 a , -a Vb 9 b , p')] dp 1 . (17) 

In practice 9 a and 9 b are taken to be 9 a and 9b from (13). The relationship between the input 
correlation coefficient and the digital covariance is plotted in Figure 5 for a fixed threshold level 
9 a = 0.61. 

For a given f a b, the correlation estimate p is determined by nonlinear inversion of (17). 
The inversion technique must be carefully chosen so that systematic errors arising from the ap- 
proximation are not larger than the statistical uncertainty of the estimate. This requirement is quite 
stringent. For example, from (5), a radiometer with a system temperature of T sys = 500 K and a 
noise requirement of A T rms — 0.1 K for the third or fourth Stokes parameter requires a measure- 
ment of p with absolute error less than 0.1K/(2 • 500K) = 1 x 10 -4 . The two existing inversion 
techniques for three-level correlators are based upon power series inversions of either the bivari- 
ate normal integral [19] or the one dimensional integral (17) [20]. In the former method [19] the 
inversion was derived for the cross-correlator, while for the latter method it was derived for the 
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auto-correlator. Both share similar convergence characteristics, e.g., third-order expansions are 
required to obtain 0.1% accuracy or an absolute error of 10 -4 for \p\ < 0.6. The latter technique, 
however, is mathematically simpler and permits an analysis of the effects of system nonidealities 
(considered in section A). Since this expression was originally derived for the autocorrelator, a 
new and more accurate expression tailored to the cross-correlator is presented here (the derivation 
is presented in Appendix A). First the integrand of (17) is approximated by a Taylor series about 
p' = 0. Next, the series is integrated to obtain: 


2 

fab — ~ exp 

n 



x 


(o, 2 -i)p 3 + 

(3 - 6 el + %) (3 - 661 + *») *> 5 ] + 0(p 7 ). (18) 


Finally, a fifth-order power series inversion [21, (3.6.25)] is carried out on (18). The resulting 
estimate maps f ab into p with absolute error of ~ 10 -5 for \p\ < 0.5 and normalized threshold 
levels 9 a , 9 b of 0.61 ± 10%: 


where 


1 - c 3 

P - —Tab ~ -4 J 
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(19) 


Cl ~ ~ exp --{91 + el) 


c 3 = =- exp 
371 


( 20 ) 


-1 W + », 2 )] (61 - 1 ) (61 - 1 ) 

C5 = 6^ 6XP [” \ ^ + *»)] ( 3 - 66 ‘ + ( 3 - Ml + »1) 

Acceptable inversion errors for Earth-science polarimetry are attainable using this fifth-order power 
series. 


B Sensitivity 

A radiometer’s fundamental sensitivity is limited by the available bandwidth, observation time, and 
receiver noise temperature. The radiometric sensitivity of a polarization correlating radiometer is: 


ATj7, rm3 = 


dp/dTu : 


( 21 ) 
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where op is the standard deviation of the estimate p. For continuous (analog) correlation using N 
independent samples and small values of p it can be shown that lim^o^p = 1 /vN [11], Thus, 
using (5) and (21) the fundamental sensitivity for analog correlation is: 

‘^‘\jT VtS y S Th t$ y S 




VN 


( 22 ) 


In the case of a digital correlator, the quantization noise of the A/D converter increases AT Urms 
above this value, the amount of which is a function of both the number of A/D converter levels and 
the threshold voltages. Using a digital correlator we have: 


AT(/ )rms — 


r ab 


drab/dTu 


(23) 


For the three level system with balanced channels (9 a = 0(, — 9), the sensitivity for vanishingly 
small correlation is (see Appendix A): 


lira A Tu,rms = 2tt [1 - $(0)] e' 

p>0 


<P y/T x 


v } ays 1 hysys 


y/N 


(24) 


The impact of quantization noise can be minimized by proper selection of the threshold voltages 
vth a ■ The optimal value of 9 (determined numerically) is 0.61 with a corresponding sensitivity of: 

A T„. rro , = 2A7 '/ T -^-«- (25) 


Comparing this expression to (22) we find that the 1.6-bit digital correlator achieves 81% of the 
sensitivity available from an ideal analog correlator. 

The total power channels are useful for normalized threshold level estimation as well as 
measurements of the first two Stokes parameters. The sensitivity of a total power channel can be 
calculated in a similar fashion by 


AT„ 


<7o2 


d(^l)/dT a 


(26) 


With the threshold levels 9 a — 9^ — 0.61 (i.e., set for optimal cross-correlator sensitivity), the total 
power channels have a fundamental sensitivity of (see appendix B) 


AT, 


a y rms 


= 2.20^ sva,a 


y/N ' 


(27) 
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The ideal (analog) total power radiometer has a sensitivity of T sys /\/N. Thus, a three-level digital 
total power radiometer can achieve 41% of the sensitivity of an ideal analog radiometer when 
the threshold voltages are optimized for the cross-correlation channel. It is noted that in (27) 
the optimal sensitivity for the total power channels is not used because the threshold voltages 
were chosen to optimize the cross-correlation channel. In otherwords, the threshold level value of 
0.61 is the optimum value for small cross correlations; however, this value is not optimal for the 
total power channels. This choice is acceptable, however, because in the polarization correlating 
radiometer the total power channels are primarily used to measure the relative threshold level 
values. If the thresholds were to be set for optimum total-power sensitivity, the digital total-power 
radiometer could 78% of the sensitivity of the analog radiometer with 9 a = 1.58. 

Ill Systematic Errors 

Two sources of systematic errors in a polarimetric radiometer are analyzed in this section: ( 1 ) errors 
caused by threshold asymmetries in the A/D converters (sampler offsets) and (2) other correlator 
gain attenuating sources: A/D hysteresis, timing skew, and phase differences. Design suggestions 
to mitigate each error are discussed. 

A Sampler Offsets 

Threshold level asymmetries in the correlator A/D converters produce systematic errors in the 
variance and correlation measurements. When extreme accuracy is not required, the effects of 
these relatively small biases can be neglected. For high-accuracy applications, such as found in 
microwave polarimetry for wind vector measurement, however, asymmetric threshold levels cause 
attenuation and offset variations that require compensation. We present here an analysis illustrating 
the second-order behavior of bias effects along with a simple analytical correction which can be 
incorporated into an in-situ radiometer calibration algorithm. The three-level A/D converter with 
a transfer function given by (8) typically exhibits small DC offset voltages Vg a , which effectively 
cause the threshold levels to be asymmetric about ground potential. We define the normalized 
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threshold asymmetry to be 6 a - v 6a /o Va . Incorporating this offset into the transfer function we 
now have (see Figure 6): 


I 1 if v > v tha 

-1 if v < —v t h a (28) 

0 otherwise 

Relations (9) and (10) can now be recomputed to reveal the effects of threshold offsets. 

A.l Correlation channel 

The digital correlation coefficient (17) including offsets becomes: 

r ° 6 = r ah \ p=0 ^a Va o Vb [f (o Va (0 a + 6 a ),a Vb (0 b + 6 b ),p') + f{o Va (0 a + 6 a ),a Vb {-0 b + 6 b ),p')+ 
f(°v a (-9a + S a ),a Vb {6 b + S b ),p') + f(a Va (-6 a + S a ),a Vh (-0 b + 5i,),p')] dp'. (29) 

Equation 29 can be considered equivalent to (17) but with small gain and offset perturbations of 
order 6 a and S b . We show here that the gain error is negligible if the input correlation coefficient is 
small. In contrast, the offset error is found to be an order of magnitude larger than the gain error. 
This correlation offset, however, is parameterized in terms of the threshold level offset and may be 
compensated via calibration using two unpolarized standards. 

The correlator offset error arises from the constant of integration r ab \ p _ 0 in (29). This 
constant was not explicitly shown in ( 1 7) because ideally it is zero. The constant can be evaluated 

by taking the expected value of (10) with p = 0 and using the modified definition of h(v — v ba ) in 
(28): 

r a6| p=0 = ( h ( v a ~ Vg a )h(v b - ^J)| p=0 (30) 

Clearly, when either threshold level is ideal (i.e., <5 Q = 0) the above term vanishes. A shift in both 
threshold levels, however, causes the offset error to become non-zero, the expected value of which 
can be separated into a product of two expected values since v a and v b are statistically independent 
when p = 0. The resulting correlation offset is: 

r oi>\p=Q = [1 ~ + 6 a ) - $(0 a + 5 0 )] [1 - ®(-0 b + s b ) - $(0 b + <J t )] (31) 
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Assuming S a and S b are small allows (3 1) to be approximated using Taylor series expansions about 
±9 a and ±0 b . The first term in the product is: 



(32) 


where the S 2 terms cancel to leave an odd-valued function. The linear behavior of (32) makes the 
threshold asymmetry a significant source of error. The constant of integration in (29) is the product 
of two such terms: 


r «»uo = +o(« 

1 1 

0 
+ 

-let 
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i i 

(33) 

= \e t ex P (Bl + dl) 

ns + 0(d a 6lS b 6 3 a ) 

(34) 

where 



3 

ll> 

<3*1 On 

b 

b 

(35) 

is the normalized threshold offset product. The threshold asymmetries thus affect the digital corre- 
lation offset by an amount proportional to the normalized offset product. Expressed using voltages: 

vs«v Sb 

ns = 

Vtha ^thf, 

(36) 


The above product is generally a slowly time varying hardware constant, but as will be shown in 
section IV, it can be estimated using a conventional two-look unpolarized calibration. 

The correlator gain perturbation is found by expanding the integrand of (29) in a three- 
dimensional power series in p\ and S b , then integrating the resulting expansion with respect to 
p'. The algebra involved (see Appendix C) is cumbersome, although the result can be expressed as 
a sum of two series. The first series r ab\g a= g b= o is the ideal relationship between p and r given by 
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(17). The second series is an error series 5r ab (5 a , 8 b , p) caused by nonzero threshold offsets S a and 
8 b . Collecting these terms we have: 

r "» = r »*Uo + + “ ex P ~ 2 + 9 *)1 x 

[(i - el) si + (i - ef) si] P + if.* (i - 2&1) (i - 2 el) <37) 

-g[K + el + 3) (i - el) <s 2 + (i - el) (ml + e; + 3) + o (p'.s 3 ) 

The above series is truncated at 0(p 4 ) and 0(8 3 ). Assuming that the nominal threshold levels are 
equal to the optimal value ( 6 a = 0.61) the error series becomes: 

8rai,(8a, S b , p) » -0.3140 (<5 2 + 8 2 b ) p + 0.0164 S a 8 b p 2 + -0.5621 (<5 2 + <5 6 2 ) p 3 (38) 

The error series is a sum of components that are 0(8 2 p), 0(8 2 p 2 ), and 0(8 2 p 3 ), respectively. 

To determine which components of the series are significant we assume that p = 0.1 and 
8 a = 6p = 8 for which the absolute values of the three components are 

|-0.3140 (8 2 a + 8 l) p\ « 0.0445 2 
|0.01645 a ^p 2 | » 0.00016J 2 
|-0.5621 (£ 2 + <5 2 ) p 3 1 » 0.00079* 2 

To render these error terms insignificant, the quantity 8 2 must be sufficiently small. Using the 
previous criterion that all errors in p of magnitude < 10 -5 are negligible, the normalized threshold 
offsets should be no larger than 1 0~ 2 , that is, vg a < 10~ 2 a Va . This is readily attainable using 
precision electronics for o Va ~ 0.5 V. If threshold offsets are not small enough, then the offsets 
should at least be controlled to render insignificant the higher-order terms (e.g., p 2 , p 3 . . .), in 
which case only a correctable gain error occurs. For this latter case, it is sufficient that v ba < 
10- V to cause the magnitude of the p 2 and p 3 terms to be less than 10 -5 . The remaining error 
is linear in p and can be modeled as an effective change in the correlator gain: 

Tab = r ab\ p= o + \ exp ~ ( 9 2 a + 9 2 h ) (l - ^ [(l - 0 2 ) 8 2 a + (l - 0 2 ) S 2 b ] ^ p + 0(p 2 , 8 3 b ) 

(39) 
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Note that the ideal correlator output ^ 6=0 is implicitly included in (39). Typically the thresh- 
old offsets are small enough so that the gain perturbation is a only few percent. 

A.2 Total power channel 

The effect of threshold asymmetry on the total power channels is a perturbed system gain and offset 
along with a residual nonlinearity that we show to be negligible. Consider the expected value of 
the total power output: 

(si) = $ (-0a + 6 a ) + 1 - $ (9 a + S a ) (40) 

This expression is a simple extension of (1 1) but includes the threshold asymmetry. The above can 
be approximated in <5 Q as 

(s 2 a ) = <* (~9*) + 1 - * ( 0 Q ) + \\j\e~ el/2 0X + O (£) ( 41 ) 

Similar to the correlation channel, the expected value of the total power channel is a sum of the 
ideal output and an additional error series. We can now show that part of the error series can be 
combined with the ideal output to compute a modified system gain and offset, with the residual 
(nonlinear) component being insignificant. If all functions of 9 a are approximated by a power 
series expansion 

$ (_0 O ) + 1 - <f> (9 a ) = 1 _ yjl 9 a + yi & l + O (91) (42) 

e" <S/2 = l -^ a + 0(0 (43) 


then (40) can be written 
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There is a gain term affecting the total power channel output by a factor of (l - l -5 2 Q ) and an offset 
of approximately This additional system gain and offset is easily identified via a standard 
two-look calibration. The nonlinear residual is 

-|Vf^ + °W)- 05) 

Assuming the optimal value for the threshold levels ( 9 a = 0.61), the above residual is found to be 
~ 10 8 a . If 8 a ^ 10 , then the nonlinear residual term becomes ~ 10 — which is insignificant 
for either total power or threshold estimation. 

B Other correlator gain attenuating sources 

Analog-to-digital converter hysteresis acts to reduce the correlation output by an amount propor- 
tional to the magnitude of the hysteresis. This effect has been modeled by D’Addario, et al. [19] 
assuming a uniformly distributed region of uncertainty about the nominal threshold. However, 
this statistical model underestimates the attenuation effect because the hysteresis is treated as a 
process that is statistically independent from the signal. Rather hysteresis is a nonstationary pro- 
cess in which the current threshold level is dependent upon the previous value of the input signal. 
To make a more accurate assessment of hysteresis a Monte-Carlo simulator was constructed to 
demonstrate the effect on the gain of the correlation channel. The simulator is based upon an A/D 
converter transfer function of the form: 

' 1 if v{nT) > v tha + v hyS!a AND h[v((n - 1)T)] # 1, 

1 if v(nT) > v tha - v hys>a AND /i[u((n - 1)T)] = 1, 

h{v(nT;v hys , a ) = ^ -1 if v{nT) < -v tha v hy3tQ AND h[v{{n - 1)T)] = - 1 , (46) 

-1 if v(nT) < ~v tha - v hy3tQ AND h[v({n - 1)T)] ^ -1, 

0 otherwise. 

v 

where Vh y3tCt is the hysteresis voltage. The transfer function is graphically illustrated in Figure 
7. Input correlation coefficients in the range —0.1 < p < 0.1 were tested with varying levels 
of hysteresis using 2 Monte-Carlo samples for each case. In Figure 8, the correlator output r 0 $ 
is plotted for values of hysteresis in the range 0 < ,s,a/&v a < 0.333. The computed relative 
attenuation on due to the hysteresis is practically independent of p over this range. The results 
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of the statistical model of [19] for the same conditions are also plotted. Compared to the Monte- 
Carlo simulations, the statistical model underestimates the attenuation effect by a factor of 10 at 

Vhys,a/ @v a = 0 . 1 . 

The results suggest that hysteresis has a moderate effect on the correlator output. A re- 
duction in correlator gain of 1 % is caused by a hysteresis voltage equal to 2% of the RMS signal 
voltage. For a 0 dBm signal into 50 this reduction would be caused by a hysteresis voltage of 
only 4.4 mV. Appropriate care should be taken in the design of the A/D converter circuitry. Alter- 
natively, in-situ correction could me made using a polarimetric calibration system or with precise 
knowledge of the hysteresis levels. 

Timing skew between the A/D converters or (equivalently) additional delay in one of the 
RF or IF paths has a similar effect of reducing the correlator output. The baseband signals at the 
correlator input can be modeled by 

v a(t) = v° a {t) + v c (t) (47) 

v b{t) = v$(t - At) +v c (t- At) (48) 

where v a (t), v®(t), and v c {t) are mutually uncorrelated and wide-sense stationary, and At is an 
additional path delay or timing skew. If v c (t) is bandlimited then the cross-correlation function is: 

Rvav b (At) = po Va a Vh sinc(2?ri? Af ) ( 49 ) 

where D is the bandwidth or bandlimiting cutoff frequency of v c (t), and the function sinc(z) = 
Thus, timing skew will reduce the measured correlation coefficient. For example, a 1% 
reduction would be caused by At = 0.039 D~ l . For a 500 MHz bandwidth, this corresponds to 
Af = 78 ps or ~ 24 mm of free-space IF path delay. 

Phase differences between the two signals will also cause an attenuation of the measured 
correlation coefficient. Phase differences can be caused by LO distribution phase errors and phase 
imbalances between amplifiers, filters, waveguides, and other transmission paths of the respective 
signals. If the phase error is denoted A 7 , then the cross-correlation at the correlator inputs is 

Rvav „ = pc r Va a Vb cos(A 7 ) (50) 
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For example, a 10° or 20° phase difference will cause a 1.5% or 6% reduction in the correlation 
coefficient, respectively. 

IV Calibration 


Calibration of a digital polarimeter entails the periodic identification of slowly time- varying system 
hardware parameters. For the total power channels, these constants are the system gain and offset. 
For the polarization correlating channel, the threshold-offset product (36) in section A. 1 as well 
as any other additive offsets (such as those originating from correlated LO noise) must be identi- 
fied. As shown below, these parameters can be estimated using the simple hot and cold views of 
unpolarized blackbody standards as obtained during conventional total power channel calibration. 


A Total power channel calibration 


Identification of the gain and offset of the total power channels in (14) allow an antenna temperature 
estimate to be made. The output of the total power channel is related to the antenna temperature 
estimate by: 


$ 


-l 



-2 


= 9a 




(51) 


where the left hand side is the linearized digital variance, g a is the radiometer system gain in K _l , 
and the receiver temperature T RE c,a is the system offset. 

The gain and offset can be estimated if the radiometer is presented two known antenna 
temperatures of differing values. The digital variance measurements corresponding to the hot and 
cold calibration temperatures, denoted and Tg l / L , are: 



-2 


-2 


— 9a ( Tqal + TfiE C ,aj (52) 

= 9a (Tcal + Trbc.o) (53) 
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( 54 ) 


This system is easily solved for g a and f RECa : 



With the two system parameters properly identified, the estimated antenna temperature as mea- 
sured by the total power radiometer is 



B Correlator channel calibration 

Estimation of the digital correlator system parameters can be achieved using the same hot and cold 
calibration looks required for the calibration of the total-power channels. There are two unknown 
calibration offsets. The first is the correlation bias po, caused by correlated LO thermal noise. The 
second is the digital correlation offset r a6 | p=0 , which is caused by threshold level asymmetry and 
is subtracted from the digital correlator output prior to conversion into the continuous correlation 
coefficient. The expected value of the correlator output given an unpolarized brightness field at the 
antenna input (i.e., T v = 0) is 

r oh Ire =0 = Cofl-* + Cl Po + Czpl 4- c 5 pl (57) 

where 

2 r i 

Co = -0«0&exp I ■ --{el + el) (58) 

Ci, C 3 , and C 5 are given by (20), and is the offset product (36). The fifth-order term c$pQ can be 
ignored ifp 0 < 0 . 1 , which is usually the case when Ty = 0 . 

The two calibration targets provide unpolarized emission at two different radiation inten- 
sities. Sequential views of the hot and cold targets provide the digital correlation measurements 
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and ?H ld for the hot and cold looks, respectively. Using these two measurements a system of 
equations can be formed: 


cfiot JhotCi Jxot'Z , hot -^3 /ro\ 

r ab = C 0 TTrfCi po + C 3 p 0 (59) 

-cxold „coldci s'Cold'Z i i&CW 

r ab ~ C 0 TTtfCi Po + C 3 Po (60) 


The coefficients cf ot and cf^ are computed by using the relative threshold values and 9^ ld , 
respectively. Using only a third-order expansion in po allows the above system to be solved ana- 


lytically. An estimate of the threshold-offset product can be found by: 




Cxold 
1 ab 


n cold C* 

C 1 Po 


n.cold'Z 3 


Po 


old 


(61) 


and an estimate of the correlation bias is a root of the following cubic: 

( Jiot \ / ,Mot \ / ,Jiot \ 

-V + ^3^“) + (<f ' - a + (<r - ^a<s") fS = 0 <«) 

The solution of a cubic equation is given in [21, (3.8.2)]. For this particular cubic there is typically 
one real root and a pair of complex conjugate roots. The real root is the desired solution for po and 
is given by: 


Po — 

where q and r are defined as 


r+{q 3 + r^Y [r - ( 9 3 + r 2 )*] * 
1 (cT - 

Q 3 (c2 ot - |£c?") 

( Shot _ Pp ot 

1 y ab c^ old T a b ) 

2 (c£ ot - ^cf d ) 


(63) 


(64) 

(65) 


Care must be taken in choosing the proper branch of the cube roots; otherwise the solution is 
straightforward. 

Once ng and p 0 have been determined, the correlation coefficient estimate p is computed as 
follows: 


P = 




(66) 
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where the estimated digital correlation bias is computed as in (A.l): 


r a b\ 


0 =o = *s B aO b exp (el + Ol) 

An estimate of the third Stokes parameter is subsequently computed using (5): 

fu = 2p\J T VjSy3 T h 


[ h,sys 


( 67 ) 


( 68 ) 


V Hardware Demonstration 


The digital correlating polarimeter concept was demonstrated using the Polarimetric Scanning 
Radiometer/Digital (PSR/D) airborne imaging instrument [2,22,23], The PSR/D is the first multi- 
channel tri-polarimetric (first three Stokes parameters) high-resolution imaging radiometer for 
aircraft-based studies of land and ocean emission (Figure 9). The radiometer operated success- 
fully in a conical scanning configuration to measure the first three Stokes parameters over the 
wind-driven ocean at 10.7 GHz (X-band) and 37.0 GHz (Ka-band) [3], In-situ calibration was 
accomplished using unpolarized hot and ambient temperature blackbody calibration targets and 
verified using a ground-based polarimetric calibration target [8]. 

A Hardware 

The 10.7 GHz radiometer was a superheterodyne single-sideband (SSB) system with a low noise 
amplifier (LNA) front-end on each channel. The 37.0 GHz radiometer was a superheterodyne 
double-sideband (DSB) system with mixer front-ends. The receiver topologies were similar to Fig- 
ure 1, with common local oscillators providing coherent downconversion of the two orthogonally- 
polarized channels. Sufficient RF and IF gain (~ 100 dB) was provided to amplify the baseband 
signals to a level appropriate for three-level A/D conversion. A separate analog total-power detec- 
tion circuit, with square law detectors and video amplifiers, was included. This redundant circuitry 
provided independent measurements of the first two Stokes parameters. 

The IF signals of the 10.7 GHz system were fed to a 1 GS/s emitter-coupled logic (ECL) 
digital correlator capable of measuring the variances and correlation coefficient of two 500 MHz 
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orthogonally-polarized baseband random signals (see Figure 10). The 37.0 GHz system provided 
a 1000 MHz wide IF band, thus, each IF band was subdivided into two 500 MHz subbands. These 
signals were fed to two correlators. By such frequency demultiplexing the entire available band- 
width of the 37.0 GHz radiometer was utilized. Each correlator comprised three functional blocks: 
the A/D converters, the multipliers, and the accumulators. The high-speed A/D converters were 
constructed using dual-window comparators yielding three levels of quantization at a 1 GS/s rate. 
The comparator threshold levels were dynamically adjusted to track slowly-varying drifts in IF 
signal power. Threshold level adjustment using eight-bit D/A converters provided (for Gaussian 
signals) a 48 dB A/D converter input dynamic range. The operating dynamic range was limited to 
~30 dB, however, by setting a minimum allowable threshold level ensured a good signal-to-noise 
ratio. The typical threshold level was ~0.3 V. 

Estimates of second-order digital signal statistics were made by squaring and cross-multiplying 
the A/D converter outputs, then accumulating using digital counters. The total power, or variance, 
of an individual channel was measured by counting the number of times the input signal exceeded 
either the positive or negative threshold levels as in (9). The correlation coefficient was similarly 
determined by separately counting the number of positive and negative correlation counts. A total 
of eight AND/NAND gates composed the entire three-level multiplier circuit. The outputs of the 
digital multiplier were accumulated in four 24-bit counters providing 16.8 ms of integration time. 
The initial 1-Gbit multiplier outputs were prescaled using high-speed 8-bit ECL ripple counters. 
The system clock was distributed differentially to the counters using 50 Q odd-mode coupled mi- 
crostrip lines. The high-speed ECL signals exhibit transition times shorter than 250 ps; therefore, 
the digital signals have spectral content >4 GHz. On-chip and interconnect propagation delays 
within the multiplier circuit were compensated with clock delays generated by programmable de- 
lay chips. To save power, the output from these counters were carried to 16-bit TTL counters. The 
most-significant 16 bits were buffered and read by computer. The circuit was fabricated on six- 
layer G10 fiberglass circuit board. Microstrip interconnects were placed on the outer two layers of 
1/2-oz copper and power was distributed via the internal layers of 2-oz copper. 
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Redundant analog total-power channels were implemented in parallel with the digital ra- 
diometers. The same IF signals fed to the digital correlators were coupled to square-law detectors 
at ~-23 dBm power level. Video amplifiers following the square-law detectors used integration 
times of 8 msec. The video amplifier output ranged from 0-10 V and was sampled by a 12-bit A/D 
converter. An analog offset was added to the video amplifier output to maintain the signal level 
within the operating voltage range of the A/D. 

B Calibration 

The unpolarized hot and ambient method of Section IV was used to identify the correlation offset 
and the threshold-offset product system parameters. Microwave foam absorber in ambient and 
liquid nitrogen conditions provided unpolarized radiation fields. The estimated system parameters 
and the radiometric sensitivities of the Ty channels are presented in Table 1 . The radiometer sensi- 
tivities were measured by computing the standard deviation of 100 calibrated measurements while 
staring at an ambient temperature absorber. The resulting measured sensitivities are 1.03-1.29 
times the fundamental limits stated in (25). The increased noise is attributed to RF and IF ampli- 
fier gain fluctuations, non-optimal selection of threshold levels due to hystersis, and threshold-level 
random noise. 

To verify the effectiveness of the unpolarized two-look technique for calibrating the Ty 
channel, a polarimetric calibration similar to that described in [8] was performed. The digitally- 
estimated third Stokes parameter Ty, as given in Section IV, can be related to an incident polarized 
field by the following: 

Tu = 2(p) \/T^^T\y^h — gwT v + guhTh + guuTu + oy (69) 

where the gains guy, gy v , and guh and offset oy are unidentified system parameters that might have 
been left uncorrected by the two-look non-polarized calibration procedure. Assuming that p 0 and 
-Ks completely calibrate the third Stokes parameter channel, then g vu = 1 and g Uv = g Uh — oy — 0 
(the gain terms are dimensionless and oy is in Kelvins). By comparing measurements of Ty with 
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an expected polarized emission vector [T v , T h , T v ] t generated by a polarized calibration target, 
the residual gain and offset terms can be checked. 

The polarized calibration target included an ambient load, a liquid nitrogen load, and a 
polarizing grid. The cold load was constructed by immersing a 56 cm x 56 cm square of the 
convoluted foam absorber in liquid nitrogen. The liquid nitrogen bath covered the absorber tips by 
at least 0.5 cm to ensure temperature uniformity. The ambient load was shrouded in a styrofoam 
jacket. The temperature difference between the targets was T hot - T cM « 210 K. The polarized 
target was mounted on a turntable that provided rotation about the feedhom axis. The rotational 
position was measured using an optical encoder with 0.25° precision. 

The details of the polarimetric calibration procedure are documented in [8], Initially, the 
radiometer antenna was aligned with the polarized target such that the incident Stokes field was 
T v = T u T h = T 2 and T v = 0. Measurements were made while slowly rotating the target 
over an angular range of ~420°. An additional piece of absorber at ambient temperature was 
used for an unpolarized look. A total of 600 radiometer samples were recorded and averaged into 
2° bins, resulting in 180 points for a full 360° rotation with ~28 ms integration time per point 
(AT R ms ~ 0.4 K). Using the calibration parameters found with the two-look unpolarized method, 
the output of the digital correlator was converted into calibrated values of f v {cti) ( i = 1 . . .n) 
for the different angles and T V (UP) for the unpolarized look. The measurement vectors and 
the columns of the observation matrices for both 10.7 and 37.0 GHz are plotted in Figure 1 1 . By 
visual inspection, T v and T h mixing into T v appears to be nonexistent, but the correlator output is 
attenuated ~25% at 37 GHz and ~70% at 10.7 GHz compared to I\j 

Using these data, the gains and offsets for the 10.7 and 37.0 GHz T v channels were cal- 
culated according to the methods in [8] and presented in Table 2. First notable are the small 
(~0. 1-0.05%) g Uv and g Uh polarization-mixing gain terms. These negligible values are a direct 
result of using the digital correlator versus a polarization adding or analog correlation design. 
Second, the offsets o v are small (-0.05 K and 0.44 K), corresponding to correlation coefficient 
offsets p 0 ~ 10 _4 -10“ 5 . Recall, this level of offset was the design goal setforth for systematic 
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errors described Section III. The negligible values of gy v , gy^, and oy demonstrate the successful 
application of the unpolarized two-look technique for calibrating the digital correlator. 

We note, however, that the correlator gains gyy are significantly less than unity. This resid- 
ual gain term cannot be identified using the unpolarized two-look technique because it does not 
provide a stimulus with known and non-zero correlation. The polarimetric calibration technique 
does provide this estimate. Sampler hysteresis and timing skew could account for ~10% of the 
attenuation, with the remainder due possibly to phase-imbalances in image-rejection and low-pass 
filters present in the radiometers. Subsequent polarimetric calibration experiments have shown that 
this attenuation does not vary by more than a few percent over several months, allowing a constant 
correlator gain coefficient to be included in the post-flight data processing routines. 

The calibration exercise also yielded analog and digital total-power measurements that 
were compared for consistency. One hundred twenty-five samples with T B = 80 to 290 K were 
compared. The mean and standard deviation of the analog-digital measurement differences are 
tabulated in Table 3. The mean of the differences /j,a-d are relatively small and range from -0.16 
to -0.62 K. The constant negative bias across all four channels suggests a slight nonlinearity of 
~0.5 K over a ~200 K range or ~0.25% full scale. The standard deviations of the analog and 
digital data are denoted a a and a D , respectively, and are equivalent to their respective AT rms ’s. 
The standard deviation of the set of differenced brightnesses oa-d is close to the root-sum-squares 
(\/ T 0 p) of the individual standard deviations. This consistency is important because it indi- 
cates that the variations in the differenced data are produced by the independant random variations 
in the individual analog and digital data sets as opposed to some correlated, external, random fix- 
ations affecting both analog and digital detectors. In all, the digital total-power radiometer tracks 
the analog system quite well. 

C Ocean Surface Emission Measurements 

Perhaps the most important application of the third Stokes parameter is the measurement of ocean 
surface wind vectors (e.g., [3]). Boundary layer winds drive the ocean surface creating gravity 
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and capillary waves and ocean foam, which in turn cause azimuthal dependencies in the upwelling 
microwave brightness temperatures [24, 25], This dependence has been observed using both satel- 
lite [26] and aircraft radiometers [27-29] in all four components of the Stokes vector. 

In March 1997, the PSR/D, along with a complement of passive and active instruments 
(together known as the Ocean Winds Imaging suite, or OWI [28]), was operated aboard the NASA 
Wallops Flight Facility P-3B aircraft as part of the Labrador Sea Deep Convection Experiment. 
Conically-scanned microwave brightness imagery of the wind driven ocean surface was collected 
during coordinated flights over the R N.Knorr. We discuss here PSR/D data collected on March 4, 
1997 from 15:00-16:00 UTC. Centered over the Knorr, the flight pattern consisted of six straight 
and level flight legs organized in three pairs, each 60° apart in heading. This “hex-cross” pattern 
covered a triangular-shaped area -30 km along a side. Data collected during the pattern consisted 
of 135 scans with 227 samples per scan at an elevation angle of 53.1° from nadir. Surface winds 
were - 14 m/s from -260° as measured by both GPS dropsonde and Knorr wind sensors. The 
ocean temperature measured at h =- 3 m was 3°C and the air temperature at h = 23 m was -12°C, 
thus, the air-sea temperature difference of 15°C indicated significant boundary layer instability. 

The measured azimuthal brightness variations at 10.7, 18.7, and 37.0 GHz are plotted in 
Figure 12. The 170 scans were averaged to accumulate a total of -2.7 seconds of integration time 
for each of the 227 azimuthal points, thus corresponding to a A T rms of -0.05 K. First- and second- 
order harmonic amplitudes were determined via a least-squares fit to the mean azimuthal signa- 
ture and are overlaid in the plot. For comparison, the 37 GHz Special Sensor Microwave/Imager 
(SSM/I) global average wind direction harmonics from [26] are also plotted with a -2 K offset 
for clarity. The plots show a distinct variation of -2-4 K with strong first and second harmonic 
dependence in the vertical and horizontal polarizations, respectively. For both polarizations, the 
dominant harmonic amplitude is larger at the higher frequency. Furthermore, the measured verti- 
cal and horizontal harmonic amplitudes at 37.0 GHz exhibit excellent agreement with the SSM/1 
global average wind direction harmonics. Although slightly lower in amplitude, the PSR 1 0.7 GHz 
wind direction harmonics are otherwise comparable to the 37 GHz SSM/I harmonics. Of particular 
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interest is the strong (~1 K amplitude) first harmonic present in the third Stokes parameter signa- 
ture. The large first harmonic content of this signature is indicative of a strong windward-leeward 
asymmetry in the ocean wave structure. The measurements are similar to simulated results ob- 
tained using an asymmetric wave geometrical optics model [24], These measurements and their 
excellent agreement with satellite and model data verify the utility, calibration technique, and per- 
formance capability of the digital polarimetric radiometer. 

VI Discussion 

The design techniques and radiometer hardware described here demonstrate the utility and tech- 
nological feasibility of the digital polarimetric radiometer for earth remote sensing applications. 
Other polarimeter topologies such as the analog correlating polarimeter or the analog adding po- 
lanmeter are possibilities. Such systems, however, can exhibit polarization cross-coupling beyond 
that caused by the antenna system that is not easily identifiable without sophisticated calibration 
techniques. On the contrary, the digital polarimeter, if built to the proper design specifications, has 
the distinct advantage of negligible Stokes parameter cross-coupling and affords in-flight periodic 
calibration of all polarimetric channel parameters. Further, use of a three-level digital correlator 

provides a simple means of calibrating both correlation offsets as well as total power measure- 
ments. 

To reiterate, the following design rules developed in Section III for the A/D converter pa- 
rameters required to limit offset and gain perturbations are: 

’• < O-Ol^ to minimize correlator offsets r ab | and Sr ab 

2. minimize hysteresis, timing skew, and phase differences to maximize g uv 

Adherence to these rules is highly desirable in order that the radiometer can be used to make 
accurate measurements of the third and fourth Stokes parameters. However, if these design speci- 
fications cannot be met, the radiometer may still be used pending regular calibration using a polari- 
metric calibration standard or a similar method such as correlated noise injection to compensate 
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for gain and offset variations. 

The PSR/D hardware demonstration confirms the ability to fabricate, operate, and calibrate 
a digital polarimetric radiometer in the field and exemplifies its utility in earth remote sensing, par- 
ticularly in the observation of ocean surface winds. As a follow-on to the PSR/D demonstration, an 
implementation suitable for satellite deployment using lower-power space -qualified CMOS logic 
at comparable sample rates is feasible and is currently under development [30], An incidental 
consequence of this work also exists in the application of digital correlators for synthetic aperture 
interferometric radiometry (e.g., [31]). Such systems will be susceptible to the same effects of 
nonlinearity, threshold asymmetry, timing skew, A/D converter hysteresis, and phase errors, all of 
which have been discussed here. 
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A Correlation Coefficient Inversion 


The digital correlation coefficient can be computed by the following: 

7*06 = ®Va®Vb I f{Pv a @ai &V),9 b ] P ) *b f{&v a @a,i &v b 9 b \ P ) ~b 

Jo 


f{-Ov a 0 a ,o Vb e b -,p') + f{-<r Va 9 a , - a Vb 6 b ; p ) dp' (70) 


Symmetry of the Gaussian probability density function allows us to write: 

fab = 2a „ 0 f f(o v J) a . o Vb 0{,', p ) + f 9 a , o Vb Oij : p ) dp 1 

Jo 

The function cr Va o Vb f(a Va 9 a , a Vb 9 b ; p) is recognized as the bivariate normal pdf [21, (26.3.1)]: 


P(0a, p) = 


" v ' 0 ’ 2n^l^ eXP [ 2(1-?) j 

Rewriting the expression for r ab using the bivariate normal pdf yields: 

fab = 2 [ p(9 a , 9 b ; p) + p(9 a , -9 b \ p') dp' (73) 

Jo 

The task at hand is to expand the integrand in a Taylor series and then integrate. The integrand of 


91 - 2p9 a 9 b + 9\ 


the above is 

I (p') = 2 p(9 a , 9 b - p') + 2 p(9 a , -9 b , p') (74) 

This can be expanded in a Taylor series in terms of p'\ 

I(p') = I(0) + I<‘> (0) p' + i /« (0) / + l/< 3 > (0) / + I/‘ 4 > (0) / + ... (75) 

The algebra involving the derivatives is quite cumbersome and the computer algebra package 
MapleV was used to evaluate the derivatives. These derivatives are 

—p(9 a ,9 b ,p) =p(9 a ,9 b ,0)9 a 9 b (76) 

°P p = 0 

P(9 a , 0 b , p) = p(9 a ,9 b , 0) ( 9\ - 1) (9\ - 1) (77) 

°f r p = 0 

SoP(9a, 9 b , p) = p(9 a , 9 b , 0) (3 9 a - 91 ) (3 9 b - 9 \ ) (78) 

°f r p = 0 

Sr P(9 a , 0 b , p) = p(9 a , 9 b , 0) (3 - 6 9\ + 9 4 a ) (3 - 6 9\ + 9 4 ) (79) 

°P p = 0 
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The derivatives of p(9 a , —9 b , p) are easily found by substituting —9 b for 9 b in the above. Because 
the first and third derivatives are odd functions of 9 a and 9 b , it is immediately seen that the Taylor 
series terms with odd powers of p will cancel leaving only the even powers of p. Adding the 
appropriate derivatives yields the following for the integrand: 

i W = MKA, o) i + \ (e\ - 1) (el - i) / 

( 3 - + ".'J ( 3 - r < + K) / + O (/) (80) 

Finally, integrating the above yields: 

r«* = \ exp [-^ (e\ + <?)] x [„+ i (el - 1 ) (el - 1 ) p’+ 

~(3~<iel+el)(3-iyel+el)p’- +o( P 1 ) (8i> 

B Digital Radiometer Sensitivity 

This appendix contains a derivation of the sensitivities of both the cross- and autocorrelating chan- 
nels of the digital polarimeter. These sensitivities are assumed to be optimized with respect to the 
A/D converter threshold level. 


A Cross-correlator Sensitivity 

The sensitivity of the third Stokes parameter cross-correlating channel is 


U.TTJXS Cl /dm 

drab/dTu 

Using the chain rule, the derivative in the denominator is expanded: 

dr ab _ dr ab dp 
dT v ~ ~di dT v 

The derivative dr ab /dp evaluated for small p can be computed using (15) and (16): 
dr 

Q- = °v a °v b [/(^A> °v b 9 b ] p) + °v b 9 b ; p) + 

f(°v a 9a, -<y Vh 9 b \ p) + f(-a Va 9 a , -a Vi 9 b \ p)\ 
= lira 2o Va o Vb [f(o Va 9 a ,o Vb Q b : p) + f(-a v J a , a Vb 9 b ; p)\ 


(82) 

(83) 

(84) 
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For the typical correlator, the channels are assumed to be balanced (i.e., o Va = a Vb = a) and the 
threshold levels for both of the input A/D converters are assumed to be equal (i.e., 9 a = 9 b = 9). 
Thus, the above is simplified: 

3t 

Tp = 55 2a2 a6; o) + /(-»». p)] 

= 2<r 2 [/(< 70 , a&\ 0) r at), 0)] (85) 

= 4o 2 f(o9,o9- 0) 

The remaining derivative is obtained from (5): 

dp _ 1 

dTrr 2\/T 71 

In the numerator of (82), the standard deviation of r a ^ is found by definition: 

a r ab = yf(ri b ) - {rat,) 2 (87) 

Under the limiting case of p -» 0, the expected value of the digital covariance is zero: 


Thus, the standard deviation now is 


lim(f o6 ) = 0. 

p— ► 0 


lim Of = lim </ (f 2 .). 

p-K) ai o-(0 V ' a6/ 


Expanding the quantity (f^) results in: 


]™(rL} = lim 

P~* 0 p — ►o 


n=l m — 1 


( 88 ) 


(89) 
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\ n= 1 m=I I 
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Evaluating the limit allows the expected value within the double sum to be separated into two parts 
because v a and v b are statistically independent when p -4 0: 


1 N N 

lfm( r aft) — 2 S S (M^(nT))M» a K))} (A(«i(nr))A(t;,(mr))) 


71=1 m=l 


(91) 
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Furthermore, because the samples are independent and identically distributed, all non-zero lag 
products are zero; { h(v a (nT))h(v a (mT ))) = 0 for n ^ m. Thus, the above double summation 
becomes 


n— 1 

= ~(hHva))(hHv b )) 

= 

Finally, by combining the previous results, the radiometric sensitivity is found: 


lim ATu rms = 

' 2a 2 f(ae,ad-,0)VN 

The above expression can be written in terms of 9 by substituting in (7) and (11): 


e 2 


.. A _ 2ne‘ 

1K n ^Tu,rms — 7-^y" 

P->° y/N 


[1 - *(#)] s/T^h, 


sys 


(92) 


(93) 


(94) 


Computing the value of 6 for the minimum AT UiTms can be done using Newton’s method. The 
optimal 9 is 0.61 with 


ATi;,rrns — / „ ; 


y/N 


sys 


Comparing this to the continuous correlator: 


y/^tSysTh t sys 

B Total-power Sensitivity 

The sensitivity of the total-power channel is found similarly: 

AT _ x/((%) 2 ) - ® 2 
d(si)/ar Ar , T ^ 

Once again, the denominator is expanded using the chain rule: 

*&> d(Z) dial) 


dT ANTa d(a % a ) dT ANTa 


(95) 


(96) 


(97) 


(98) 
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The first term in the product is the differential relationship between the input voltage variance and 


the output of the total-power channel of the digital correlator. From (11), 

3«) d «) [ 

<V 2 ? P [ 2 \ a ,J 

Recalling (14), the variance of the input voltage signal is 


a v a — RokBG a (T ANT ,a + TreC,o) 


( 100 ) 


so that 


3«) _ < 

9T ANTa T 3 y 3a 

where T sy3a — TnEC,a + T A nt,cc Combining the two derivatives yields 


9<?a) _ 

9TANT,a T 3y3 


1 ( Vth,a A 

iQ y/2n V °v a ) 


-- 9 n e 2 0 “. 


exp - 


1 ( V t h,a 


2 V a„ 


T sys , a V^ Q ' 

In the numerator of (102), the standard deviation of 3^ is found by definition: 

\A(3|) 2 > - (Z) 2 

This is computed by first expanding the expected value of (s^) 2 : 

({sl) 2 ) = ( -^^^(vainT)) \ 

\ L n=l J L 771= 1 J / 


= A?2 5^ S ( h 2 {v a {nT))h 2 (v a (mT ))) 


( 101 ) 


(102) 


(103) 


(104) 


Because the samples of v a are independent and identically distributed, the above expected value 


can be written: 


(HWWMrnT))) = { ^ntZ <'» 5 > 
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Conveniently, by the definition of h(v), h 4 ( v ) = h 2 (v). Substituting these two cases into the double 
summation yields 


{(3D 2 ) = iE {h 2 {v a {nT ))) + ^2 (/i 2 MnT))) </i 2 (u Q (mT))> 


TV 2 ^ " / ' tv 2 . 

n=l n=l m=l 


Thus, the variance of is 

v^) 2 ) - <^> 2 = + (l - {^> 2 - < 5 ^ 


= (1 (®a)) 


( 106 ) 


(107) 


= A /-[i-$(0 a )]$(<y 

Forming the quotient (97) with the variance (107) and the derivative (102) produces 


AT„ 


\/27r i 


e 5 V* (^a )-^ 2 Wa) 


1 

//V ' 


(108) 


The optimum value of 0 a for best total power radiometer sensitivity is found by minimizing 
AT a rtns with respect to 0 n . Again, this can be done numerically. The sensitivity is minimized 
at 0 Q = 1.58 such that 


AT„ 


= 1.28^P^ 

y/N 


(109) 


For operation of the polarization correlating radiometer, however, the value of 9 a is set to minimize 
the noise of the T v measurement. At the optimal point for T v where 9 a =0.61, the sensitivity of 
the total power channel is 


AT, 


a,rms 


— 2,20 ^ sy3 ’ a 

y/N 


( 110 ) 
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C Threshold Offset Effects 


Perturbations in the threshold voltages will change the output of the digital correlator. With asym- 
metric threshold levels, the digital correlation coefficient can be computed by the following: 

r a b = r ab | p=0 + o Va o Vh f f(cr v J6 a + S a ), a Vb (0 b + S b ) ,p') + f (a Va (0 a + 6 a ),cr Vb (-9 b + 6 b ), p') + 

J 0 

f (<7v a (-0 a + 5a ),o Vb (0 b + fi b ) ,p') + f (a Va ( -9 a + <*„), o Vb (~9 b + 6 b ) ,p')dp'. (Ill) 

Using the definition for the bivariate normal pdf (72), the above becomes 
r ab = r ai>| p= o + f p{9 a + 6 a , 0 b + $b, p') + p{9 a + <^ 0 , —9b + fib, p') + 

J o 

p(- 0 a +fia,Ob + fib, p') + p{-9 a + fia, -0 b + fib, p')dp'. (112) 

To determine the behavior of r ab with respect to fi a and fi b , we need to express the integrand in a 
three-dimensional power series in p, fi a and fi b . The integrand is 

I (p) = p(0a + fia, @b + fib, p) + p{9a + fia, ~O b + fib, p) + 

P(~Oa + fia, 9b + fib,p') T p{ 9a + *a, ~9 b + fib, p') (1 13) 
and its expansion in p' is nearly identical to that in section A. The Taylor series expansion is 

/ (p') = / (0) + /(■) (0) (0) / + I/W (0) / + A/«> (0 )/ + ... (114) 

What remains is to determine the two-dimensional Taylor series in fi a and fi b of each of the terms of 
the integrand. Throughout it is assumed that fi a and fi b are 0(6) and the Taylor series are generally 
truncated beyond O (fi 2 ) . 

The first term of I (p') is 

I (0) — P(@a + fia, 0 b + fib, 0) + p(9 a + fia, ~9b + fib, 0) + 

P(~0a + fia, 0 b + <M) +P(~9a + £«, ~9 b + S b , 0) (115) 

There are four instances of a term with the form p(x + a, y + b, 0), where x = ±9 a> y = ±6 b , 
a = fi a ,b = 6 b . The two-dimensional Taylor series about x and y will be applied to the above. 
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The expression p(x + a, y + b, 0) is the product of two standard normal pdfs: p(x + a,y + 
6,0) = Z(x + a)Z(y + b ), where Z(x) = ^=e _l2 / 2 is the standard normal curve [21, (26.2.1)]. 
The Taylor series of this product is 

p(x + a,y + b, 0) = Z(x)Z{y ) + Z {1) (x)Z{y)a + Z(x)Z^(y)b + 

~ [. Z^{x)Z(y)a 2 + 2Z^(x)Z^(y)ab+Z(x)Z^(y)b 2 ] +... (116) 

The derivatives Z 60 (a;) of the standard normal pdf can be expressed in terms of Hermite polyno- 
mials [21,(26.2.32)]: 

Z {n) (x) = (-l)-»2 - n ' 2 Z(x)H n (^=), (117) 

where the first three Hermite polynomials are H 0 (x) = l,H 1 (x) = 2x, and H 2 (x) = 4a: 2 - 2. For 
example, the second derivative of Z(x) is 

ZW( X ) = (_l)- 2 2 ~ 2 ' 2 Z(x)H 2 (-^=) 

v 2 (118) 

= Z(x) (x 2 - l) 

Substituting these results into (115), the first term of the integrand is 

/ (o) = z(e a )z(e b ) + zW(e a )z(9 b )5 a + z(9 a )z^(9 b )s b + 

i [^ (2 H0«)^(^)^ 2 + 2Z( 1 )(0jZ( 1 )(0 6 )<S a 4 + Z(0 a )Z( 2 )(0 6 )ci 2 ] + 

Z(9 a )Z(9 b ) + Z^(9 a )Z(9 b )S a - Z(9 a )Z^(9 b )S b + 

~ [zW(9 a )Z(9 b )6l - 2 Z^(9 a )Z^(9 b )S a 6 b + Z(9 a )Z^(9 b )6 2 } + 

Z(9 a )Z(9 b ) - Z^{9 a )Z(9 b )6 a + Z(9 a )Z^{9 b )8 b + (J 19) 

| [zW(9 a )Z(9 b )6l - 2 Z^(9 a )Z^(9 b )6 a 5 b + Z(9 a )Z™(9 b )6 2 ] + 

Z(9 a )Z(9 b ) - Z^(9 a )Z(9 b )S a - Z(9 a )Z^(9 b )S b + 

% [Z {2) (9 a )Z(9 b )6 2 a + 2 Z^(9 a )Z^(9 b )6J b + Z(9 a )Z^(9 b )S 2 } +... 

= 4 P(x, y, o) (l - ^ [(l - 9l) 6l + (1 — 9 2 ) tf 2 ] j + 0(S 3 ) 
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The second term in the Taylor series of I ( p ') is 


l{1) (U)p'= ^-p(0 a + S a ,e b + 6 b ,p) I + ^-p{9 a + 8 a ,-9 b + S b ,p) + 

d Q | 

Q-p(-0a + S a ,0b + S b ,p) + — p{-9 a + 6 a , -9 b + 6 b , p)\ p' (120) 

p=0 P | p — o 

The first partial derivative with respect to p , evaluated at zero, of the bivariate normal pdf is 


-p(9 a , 9b, p) = p{9 a , e b , o )9 a e b 

p=0 

= 9 a Z{9 a )9 b Z{9 b ) 


( 121 ) 


Thus, 


/ (1) (0) = (0 a + 6 a )Z(9 a + 6 a )(9 b + 5 b )Z[9 b + 6 b ) + 

{9 a + S a )Z(9 a + 6 a )(-9 b + 8 b )Z(-0 b + 6 b ) + 

( 9a 4- 8 a )Z{—0 a + 5 a )(9 b + 8 b )Z(0 b + 8 b ) + 

( 122 ) 

{-9 a + 5 a )Z(-6 a + 6 a )(-9 b + 6 b )Z(—9 b + 8 b ) 

= [{9 a + 8 a )Z(9 a + S a ) + {-9 a + 6 a )Z(-9 a + 5 0 )] x 
[{9b + S b )Z(9 b + 5 6 ) + (~9 b + 8 b )Z{-9 b + ^)] 

Next, consider the expression ( 9 a + 5 a )Z(0 Q + 8 a ) + {-9 a + 8 a )Z(-9 a + 8 a ) by expanding the 
standard normal pdfs in one dimensional Taylor series: 


{9 a + S a ) Z(9 a ) + ZU(9 a )8 a + ±Z( 2 \9 Q )6 a + ...] + 

(- 9 a + s a ) Z(-0 Q ) + Z^(-0 Q )S a + iz^(-9 a )8 2 + ... 

= 9 <* [ Z ^ + z{1] (9 a )6 a + ^ ZW{9 a )8l + ... - 9 a Z(0 Q ) - Z^(9 Q )8 a + ^Z^(9 a )8 2 Q + ... 

+ j Z(9 a ) + zW(0 a )6 a + lzW(0 a )8 2 Q + ... + 8 a Z(0 Q ) - Z^{9 a )6 a + ~Z^(9 a )8 2 + ... 

— 9 a (2 Z^ \ 9 a )5 a + ...) + 28 a Z(9 a ) + y Z^(9 a )8 2 + ... 

= Z(9 a )6 a (1 - 2 9l) + 0 (<£) 

(123) 
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Finally, substituting the above series into the a and b terms of (122) yields the following: 


/ (1) (0) = 2 6 a 6 b p(9 a , 9 b , 0) (1 - 2 0 2 ) (1 - 2 9l) + 0(6 3 ) (124) 

The third term in the series expansion of / ( p ?) contains 
/ (2) (0) = [(6 a + S a ) 2 - 1] Z(0 a + 6 a ) [(0 b + 6 b ) 2 - 1] Z(9 b + <S 6 ) + 

[(@a + 6 a ) 2 — l] Z(9 a + S a ) [(— 9 b + 6 b ) 2 — l] Z(—9 b + 6 b ) + 

[(-0a + 6 a ) 2 - 1] Z(-9 a + 6 a ) [(9 b + 6 b ) 2 - 1] Z(9 b + 6 b ) + 

r (125) 

[M« + s a ) 2 - 1] Z(-9a + 6 a ) [(-9 b + 6 b ) 2 - 1] Z(-9 b + 6 b ) 

= [(*« + 2 0 a 6a + 6 2 a - 1) Z(9 a + 6a) + (9 2 a - 2 9 a 6 a + 61 - 1) Z(-9 a + <5 a )] X 

[(9 2 + 2 9 b 6 b + 6 2 - 1) Z(9 b + 8 b ) + (9 2 - 2 9 b 6 b + 6 2 - l) Z(-9 b + 5 b )} 

Consider the following expression by expanding Z in a power series: 

(9 2 a + 2 9J a + 5 2 a - 1) Z(9 a + 5 a ) + (9 2 a - 2 9J a + S 2 a - l) Z(-9 a + <5 Q ) 

= (0l + 29 a 6 Q + 6 2 Q -l) Z(9 a ) + Z^(9 a )6 a + ~Z^(9 a )6 2 a + ... + 

(9 2 a - 2 0J a + 61-1) Z(-9 a ) + ZV(-9 a )6a + ^zW(-9 a )6 2 a + ... 

L , ' J (126) 
= 2 K + % - 1) [Z(0a) + yZ( 2 \9 a )6l + ...j + 4 9J a [Z^(9 a )6 a + ...] + 

= 2 (0l + tl- 1) Z(9 a ) + ±Z(9 a ) (el - 1) 61 +8 9 a 6 a e a Z(9 a )6 a + ... 
^Z(9 a )((9l-\)6l+(ml + 9l + Z)6l + 29l-2) + ... 

Substituting the above result into (125) produces, 

/ (2) ( 0 ) = p (0 a , e b , 0 ) [(60 2 + e A a + 3) <5 2 + 20 2 - 2 ] [(60J 2 + 6^ + 3) 6 2 b + 29 \ - 2 ] + O (<5 4 ) 

(127) 

Combining and integrating the first three terms of the integrand (119), (124), and (127) 
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results in 


r ab — r ab\ p =Q + 4 p{O a , ft> 0) ^1 - [(l — 0 2 ) 5% + (l — 0j) p + 

ftftp(ft, Ob, 0) (1 - 20 2 ) (l - 20j) p 2 + 

(0a, ft, 0) ((60* + 0 4 + 3) J 2 + 20 2 - 2) ((60 2 + 0 4 + 3) 6 2 b + 2 0 2 - 2) p 3 + 

0(p 4 ,S 3 ) 

= r ab \ p = o + 4p(0 a , Ob, 0) P + ^ (0 a 2 - 1) (ft 2 - 1) p 3 ] + 

4pW» ft, 0) [-1 [(1 - 0 2 ) 5 2 + (1 - 0 2 ) Jj] p+ 

-S a db (l — 20 2 ) (l — 20 2 ) p 2 + 

- J IW + d * + 3) (1 - ft 2 ) ft 2 + (1 - 0 2 ) (60 2 + 0i + 3) 6*] p 3 j + 

0{p 4 ,6 3 ) 

(128) 

There are two different power series in p that can be identified in the above. These are the ideal 
relationship between p and r and an error series caused by nonzero threshold offsets 5 a and S b : 

2 r i i 

rab = r ab \ p=0 + r ab \ Sa=Sb=Q + - exp -- (0 2 + 0 2 ) x 

[(1 - «l) #. + (1 - ol) dj] p+ idA (1 - 2 el) (1 - 20*) ?+ (129) 

-^1 + *1 + 3) (1 - K) £ + (1 - «?) (Wj + ol + 3) ,S*]„ 3 ] + O (p\6 3 ) 
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Table 1 : Digital correlating polarimeter system parameters found using hot and ambient unpolar- 
ized calibration targets. 


Table 2. Gain and offset terms for the 10.7 and 37.0 GHz digital polarimeters as measured using 
the polarized calibration standard. 

Table 3. Comparison of digital and analog total-power radiometer measurements. 

Figure 1: Block diagram of a typical digital polarimetric radiometer. This direct correlating po- 
larimeter, utilizes a dual polarized antenna, dual channel superheterodyne receiver, and a 3-level 
digital correlator. The IF signals are also coupled to traditional square law detectors and video 
amplifiers. 


Figure 2: Ideal transfer function of three-level A/D converter. 

Figure 3: Digital logic for the multiplier and squaring operations. The values of the boolean input 
variables A + , A~, B + , and B~ are the outputs of the 3-level A/D converters (see above table). The 
outputs C A and C B are the count enable signals for the total-power accumulators and C + and C_ 
are the count enable signals for the positive- and negative-correlation accumulators. 

Figure 4: The digital variance as a function of input RMS voltage at a fixed threshold level. The 
optimal threshold level for correlation 6 a = 0.61 occurs at a Va /v tha = 1.64 (•) and for total power 
measurement (+) d a = 1 .58 occurs at o Va /v tha = 0.633. 


Figure 5: The digital covariance r ab versus the input correlation coefficient p for 6 = 0.30, 0.61, 
1.0, and 1.58. 

Figure 6: Transfer function of three-level A/D converter with threshold offset v s . 

Figure 7: Transfer function of three-level A/D converter with hysteresis magnitude Vh ya - 


Figure 8. The reduction in the digital correlator output as a function of hysteresis amplitude. The 
results of the Monte-Carlo model (solid line with open circles) shows that the statistical model [ 1 8] 
(solid line) significantly under estimates the effect. 

Figure 9: (a) Polarimetric Scanning Radiometer (PSR) installed in the bomb-bay of the NASA Wal- 
lops P-3B Orion (N426A). (b) PSR scanhead with antenna lenses visible. The dual-band 10.7/37.0 
GHz antenna is indicated. 

Figure 10: Digital correlator module for Polarimetric Scanning Radiometer (PSR/D). The correla- 
tor is constructed using ECL discrete logic on a six-layer microwave circuit board. This module 
contains two complete correlators (including total power channels) on one board. 
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Figure 11: Plots of (a) 37.0 GHz and (b) 10.7 GHz polarized target measurements Ty(a ) (denoted 
by crosses) and the observation matrix. The columns of the observation matrix are plotted as 
follows: T v (a) (dash-dot), T h (a) (dotted), and T v {a) (solid). 


Figure 12. PSR azimuthal harmonics exhibiting wind direction dependence of the first three Stokes 
parameters at 10.7, 18.7, and 37.0 GHz. Data for Ty at 18.7 GHz was unavailable. The wind speed 
was 10 m s" 1 ±0.6 m s" 1 from 1500 to 1600 UTC on March 9, 1997. The solid lines represent the 
reconstructed second-order harmonic expansions and the dashed lines are the ±lcr r error curves 
for 170 scans. Individual points indicate mean measured brightness deviations. The 37.0 GHz 
SSM/I global average wind direction harmonics, denoted by solid lines, are shifted by -2 K for 
clarity. 
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Parameter 

X-band (10.6-10.8 GHz) 

Ka-band (36-38 GHz) 

\J Trec , vT REC ,h (K) 

1394 

939 

Pi o 

0.020 

0.078 


0.0024 

0.0058 

Receiver type 

LNA/SSB 

mixer/DSB, two subbands 

ATt/^iTheory (K) 

0.51 

0.35 

AT)/ jrmj ,:M eas’d (K) 

0.66 

0.36 


, 



- 
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